Spin Hall effect of conserved current: Conditions for a nonzero spin Hall current 
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We study the spin Hall effect taking into account the impurity scattering effect as general as 
possible with the focus on the definition of the spin current. The conserved bulk spin current (Shi 
et al. [Phys. Rev. Lett. 96, 076604 (2006)] ) satisfying the continuity equation of spin is considered 
in addition to the conventional one defined by the symmetric product of the spin and velocity 
operators. Conditions for non-zero spin Hall current are clarified. In particular, it is found that (i) 
the spin Hall current is non-zero in the Rashba model with a finite-range impurity potential, and 
(ii) the spin Hall current vanishes in the cubic Rashba model with a 5-function impurity potential. 
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Spintronics is one of the most promising new technolo- 
gies, where the spin degrees of freedom of electrons in 
semiconductors are manipulated and utilized for func- 
tions such as memory, operation, and communication. 1,2 
One of the key routes to spintronics is to invent an ef- 
ficient method to inject spins into semiconductors. In 
this respect, the intrinsic spin Hall effect (SHE) has at- 
tracted recent intensive attention since its theoretical 
proposal. 3 ' 4 This can give a larger effect by orders of 
magnitude than the extrinsic one based on the impurity 
scatterings proposed long beforei 5 i 6 i 7 

Recently two experiments have been reported on ob- 
servations of the SHE in GaAs and related materials 
Kato et al& observed the Kerr rotation due to the spin 
accumulation (~ 10/i_B(£im) -3 ) near the edges of the n- 
type GaAs sample. They suggested the extrinsic mech- 
anism of the SHE since it was almost insensitive to the 
crystal orientation. Wunderlich et al& observed the cir- 
cularly polarized LED signal from spin-polarized interfa- 
cial two-dimensional holes in p-type GaAs system. From 
an estimation of transport lifetime, they concluded that 
the spin accumulation is due to the intrinsic SHE. 

The debates on the impurity effect on the intrinsic 
SHE have continued, which are in parallel to those for 
the anomalous Hall effect in ferromagnets. In the latter 
case, the intrinsic mechanism of Karplus-Luttinger— was 
criticized and extrinsic mechanisms due to impurity scat- 
terings were proposed^ The Hall conductivity is a sin- 
gular function as the disorder strength approaches zero. 
In the metallic case, the vertex correction in the diagram- 
matic language incorporates a deviation of the electronic 
distribution function from equilibrium, and it represents 
the dissipative current. This situation is similar also in 
the spin Hall current. 

Actually, the disorder effect on the spin Hall current 
of the Rashba model in two dimensions has been in- 
tensively studied. Sinova et al£ obtained the universal 
value e/(87r)i 2 i for the spin Hall conductivity (SHC) <r s H 
without disorder. When the self-energy correction due 
to impurity scattering is taken into account, the SHC 



afj is reduced continuously as a function of the disor- 
der strength from the universal value e/(8ir)Mi On the 
other hand, Inoue et alm^ studied the vertex correction 
and found that a s H vanishes in the clean limit within the 
Born approximation. Furthermore, it has been shown 
that the SHC a s H vanishes for any value of the lifetime 
r using the Keldysh formalism^i^iii the Kubo formula 
analyticallyfiSii2iS£ and numerically, and the Boltz- 
mann equation.- 23 Thus after long debates, people have 
reached the consensus that the SHC a s H for the Rashba 
model vanishes for any 

However, the vanishing result of o~ s H depends on the 
definition of the spin current. In the previous calcu- 
lations on the SHC, the spin current is defined ad hoc 
as a symmetrized product of the spin and the velocity 
J s = \{v y ,s z }i where v y — dH/dp y , in response to the 
electric field E along the x axis. However, this "conven- 
tional" definition of spin current loses its physical founda- 
tion when the spin-orbit coupling is present and the con- 
servation law of the spin is violated. This is an important 
issue since the concept of "current" depends crucially on 
conservation; non-local effects of the current comes from 
the fact that an incoming flow goes out without loss. 
Therefore we need to search for a proper definition of 
a conserved spin current in the bulk. From this view- 
point, the conserved spin current deserves scrutiny. 
If J dV(s z ) = 0, it satisfies the Onsager's reciprocity rela- 
tion and the continuity equation for the spin in the bulk: 
dt{s z ) J r'^-Js = 0. Furthermore it is free from an artifact 
that the spin current is proportional to time derivative of 
the spin operator ji&iS hence, the spin Hall current can 
be nonzero even for the Rashba model. Even though the 
experiments on the SHE up to now^ detect the spin 
accumulation at the sample edges, we stick here to the 
SHE defined by the bulk spin current. This is because we 
consider that the generation of the conserved spin cur- 
rent in the bulk is a more fundamental phenomenon than 
the spin accumulation at edges, which is not solely deter- 
mined by this spin current since the continuity equation 
for the spin is not satisfied there. In principle, there 
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should be other means to detect the SHE without using 
the spin accumulation such as voltage measurement with 
the injected spin current^ 

In this paper, we study the SHE as generally as possi- 
ble taking into account the disorder with the definitions 
of the conventional and conserved spin current. Apply- 
ing this consideration, some new results are obtained for 
the Rashba and cubic Rashba models. We employ the 
Keldysh formalism by which the infinite series of 
the Feynman diagrams both for the self-energy and ver- 
tex correction are taken into account compactly, and the 
expression for the SHC is obtained for both definitions of 
the spin current described above. 

We consider a generic model with spin-orbit coupling 
with a random impurity potential, this random potential 
is assumed to be spin- independent, and the time-reversal 
symmetric model exclusively. In the Keldysh formalism, 
a Green's function matrix G is introduced, 
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(b) 
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where the superscripts R, A, and < denote the retarded, 
advanced and lesser Green's functions, respectively. The 
self-energy matrix £ is defined similarly. The Green's 
functions satisfy 

(G 11 )- 1 <g> G< - G< ® (G- 4 )- 1 = X<®G A -G R 

(2) 

(Go 1 - S fi » A ) <g> G R ' A = 5(1-2), (3) 

where {A®B)(l, 2) = J d3j4(l, 3).B(3, 2) and Go is the un- 
perturbed Green's function. We then separate the center- 
of-mass and the relative coordinates and perform the 
Fourier transform to the relative coordinatesi^LSi The 
final result is written in terms of the center-of-mass co- 
ordinates (T,R) and the relative momentum (w,p). We 
put the constant electric field E — (E, 0, 0) and look for 
solutions independent of T and R. Therefor the quantum 
Boltzmann equation (QBE) is written by 

[G, H] = -icE • V P G - ^eE • {V p H, O^G} 

-~eE • ({V P S, d u G} - {&,£, V P G}) + [£, G].(4) 

Here one can show from Eq. Q that the time deriva- 
tive O of an arbitrary operator O, which is independent 
of p, R, T and uj, has a vanishing expectation value in 
the steady state, (O) = 0, even with a general form of 
the impurity potential>i2iiS We start with 
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I J 2lT J (27T) Z 



du) f d 2 p 



2tt J (2?r) ; 



■tr(C[G, H] 



(5) 



where [ ] < denotes the lesser (upper right) component 
of the matrix in the Keldysh space. Eq. (JIJ is plugged 
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FIG. 1: Diagrammatic representation of the self-energy E in 
the present self-consistent approximation. Doubled line is the 
dressed Keldysh Green's function G_. Dashed line with a cross 
denotes the average over the impurity positions in the second 
order for (a) and in the third order (b), respectively. 



into Eq. (J5J, and evaluate the respective terms in the 
RHS of Eq. J3J. The first term zeE ■ V P G < becomes zero 
after an integration over p. The second term vanishes 
after an w-integration. The third and fourth terms van- 
ish after partial integrations in terms of p and u). Lastly, 
to evaluate the last term ([S, G] < ), we need a relation- 
ship between the self-energies and Green's functions. We 
employ the self-consistent Born approximation [the dia- 
grams in Fig.^(a) and (b))SS for the impurity scattering. 
Up to the second order it is given by 

72 „/ 



S(a;,p)=n< J |Vp, p /| 2 G(w, p' 

d 2 p' d 2 p" 



V^/VA P')G(w, p").(6) 



where rii is an impurity density, and Vp, p < is the Fourier 
transform of the impurity potential. From the second 
Born approximation Eq. jnj) one can easily show (O) = 
for arbitrary forms of impurity potentials? . This holds 
true even for higher-order Born approximation. On the 
other hand, for the charge current v = R, this argument 
does not apply, and (v) can be nonzero in the steady 
state as expected. 

Now we proceed to a closed form of the SHC. The 
expectation value of the conventional spin current J s = 
\{v y , s z } is obtained as 



E 



i f dio d 2 p 
hm / — - — -rtr 

e^o IE J 2tt (2tt) 2 



[dPy 



(7) 



where [ ]e implies that we retain the terms linear in a 
uniform electric field E = Ex. 

Let us turn to the second definition of the spin 
current J s as proposed by Shi et al£& J s is defined 
to satisfy dt(s z ) + V ■ J s — 0, and is divided into 
(Js) + P T , where J s = ^{v,s z }. The second term 
P T is called the torque dipole density, and is required 
to satisfy limn-^o lim^ ((s 2 (il, Q)) + iQ ■ P T (Q,Q)) = 

0, where (il, Q) are the Fourier components of the 
center-of-mass coordinates. We put Q x = 0, and 
take the limit Q y = Q — > 0, i.e., P T = P T ■ y = 
- limo^o liniQ^o ^(sz(^: Q))- More explicitly, E is spa- 
tially modulated along the y-axis, E = Ee l< 3 Y ~ tnT x 25 
Note that (s z (0, 0)) = from Eq. JSJ, which means P T is 
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finite and well-defined. In response to the electric field, 
the Green's function G and the self-energy E acquire 
terms proportional to e l Q Y - lQT , p r [ s expressed as 



1 f dw 
-i lim — / — 

Q^o Q J 2tt 



d 2 p 



tT( Sz [H 7 G}<). (8) 



We first write down the QBE to the linear order in E and 
to the linear order in 12 or Q. Next, we replace the term 
[H,G] < in Eq. (JSJ with the corresponding term in the 
QBE. While there arise a number of terms, most of them 
give no contribution to P T after partial integrations over 
Pi or to. To calculate the remaining terms, we note that 
this electric field necessarily accompanies a magnetic field 
according to the Maxwell equation. To deal with the re- 
sponse to these two fields, it is convenient to consider the 
corresponding vector potential, A = (Ae l Q Y ~ inT , 0, 0). 
Relevant terms in the QBE are classified to those propor- 
tional to itlA = E (electric field) or those proportional to 
— iQA = B (magnetic field). The resulting form is a sum 
of the contributions from the response to a dc electric 
field Ex and that to a dc magnetic field Bz: 



Js ,. i 
—— = lim — 

E e^o E 



lim — 

b^o B 



duj d 2 p 
2^ (2tt) 2 

dbj d 2 p 



tr 
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f dE 

\d Py 



G 



2tt (2tt) 2 



rtr 
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E 

^= g" 

duj ' — 



(9) 
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Here J s = J s ■ y, and [ ]g retains the terms linear in 
an uniform magnetic field B = Bz. Because B does not 
drive the system off-equilibrium, the relation in equilib- 
rium, G< = (1 - tanh(cj/2fc s T))(G^ - G^)/2, is sat- 
isfied even in the presence of B = Bz. The details of 
the derivation of Eq. will be presented elsewhere^. 
Remarkably, in calculating the total conserved spin cur- 
rent J s = (J s ) + P T , there appears a term in P T which 
exactly cancels (J s ). We also note that this formula is 
quite generic and applies to any models. The expression 
of the charge current is obtained just by replacing s z by 
— e in Eq. @, and the £?-term in the final formula is 
reminiscent of the Streda formula^SiS 

We now discuss the explicit models based on the results 
obtained above. Here we condider the Rashba model 



H 



2m 



\(cr x p) • z + v, 



and the cubic Rashba model 



H 



iX 



2m 



(p 3 _a+ -p 3 + <T-) 



(10) 



(11) 



Here, <r = (cr x ,cr v ,cr z ) = 2s is the Pauli matrix, p± = 
p x ± ip y , <j± — <j x ± ia y and v is an impurity ran- 
dom potential. We take the unit where fi = c = 1. 
The Rashba model represents an n-type semiconductor 
in two-dimensional heterostructure. The second term 
in Eq. (|10|l represents the spin-orbit coupling with an 
inversion-symmetry-breaking potential along the z direc- 
tion perpendicular to the plane. As noticed in Refs. 



and [II, J s is proportional to s y : J s = - 4imX . 

Therefore from the generic argument in Eq. (JSJ), there 
occurs no spin Hall current when this definition of the 
spin current is employed for the Rashba model (Table 
1(a)). This result is consistent with previous works us- 
ing various methods, including the calculations in the 
clean limit by Kubo formulaji^2& and Keldysh formal- 
ismiA& For finite r (cft 3> 1), it is also consistent with 
the analytioi^iiS and numericalSiiS 2 ! results by Kubo for- 
mula, and with the results by Keldysh formalismii£ii£ 
The calculation in Rcf. 19 by Kubo formula for finite r 
(e^r 3> 1) is similar to ours by the Keldysh formalism. 
Nevertheless, our approach better reveals the reason why 
the spin Hall current generally vanishes when J s oc s y . 

The cubic Rashba model, on the other hand, describes 
the heavy-hole bands of cubic semiconductors in het- 
erostructute. The conventional spin current J s can no 
longer be expressed as O; hence the previous argument 
for vanishing spin Hall current does not apply. Indeed, 
the resulting SHC is nonzero^ [Table 1(b) ], which is 
consistent with Refs. |3^ and [3^. Therefore the zero or 
nonzero spin Hall current for the conventional spin cur- 
rent J s is mostly determined whether it is expressed by 
the time derivative of some local operator or not. 

Next we turn to the conserved spin current J s . Both 
in the Rashba and the cubic Rashba models, the B-term 
in Eq. 10 vanishes, because the Hamiltonian lacks a a z - 
term, and the self-energy is independent of the spin. The 
remaining _B-term in Eq. (J5J) for the Rashba model is 
calculated as follows. In the first Born approximation 
[the first term in Eq. JHJ], J s = for general V p - P '. In 
the second Born approximation, we obtain 



Js = — 



■tr 



■G(lu,p')G(lo, P ")} 



duJ d 2 p d 2 p' d 2 p" 
2^ (2tt) 2 (2tt) 2 (2^p 

<9(WV-p"V 



dp' y 



(12) 



Even for higher-order Born approximation, the formula 
for Js can be written down. We can then see that J s 
always depends on d Py V, and the spin Hall current for 
J s is extrinsic for both the Rashba and the cubic Rashba 
models, depending explicitly on the impurity potential. 
All these considerations are summarized in Table I. With 
the (5-impurity potential the conserved spin Hall current 
vanish both for the Rashba model and the cubic Rashba 
model. We note that calculated the conserved SHC with- 
out disorder is a s H = e/87r for the Rashba model and 
ajj = — 9e/87r in the cubic Rashba models These re- 
sults are reproduced in our calculations by neglecting the 
self-energy of the lesser Green's function and taking the 
clean limit. 

We now demonstrate the finite spin Hall current for 
Js in the Rashba model with the short-range (not 5- 

function) impurity potential V(r) — Ue~^ with j3A <C 
1, where U is the magnitude of the impurity potential, f3 
is the size of the potential range and A is a momentum 
cutoff. As an approximation, we substitute the Green's 
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(a) Rashba model 



Impurity potential 


Born approx. 


Definition of spin current 
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1st 









higher 





Finite 


(b) Cubic Rashba model 


Impurity potential 


Born approx. 


Definition of spin current 






W 


Js 


«5(r) 


lst/higher 


Finite 35 ' 36 





Vp-p' 


1st 


Finite 







higher 


Finite 


Finite 



TABLE I: Spin Hall effect in the (a) (linear) Rashba model 
and (b) cubic Rashba model, with various types of the im- 
purity potential for the two different definitions: (J s ) is the 
conventional spin current, J s = ^{v y ,s z }, and J s is the con- 
served effective spin current, J s — (J 3 ) + P T . We show the 
new results in the boldface. 




FIG. 2: Contour map of the SHC o s H . It has the maximum 
value nearly along a = 2Xpft — 1. The momentum cutoff is 
set to be A « 2>Pf, and U/ef = 0.02. 



functions in Eq. with those for the <5-function impurity 
potential, and calculate J s . Figure [3 shows the results 
for the SHC <j s h in the parameter space of Xpf/sf and 
tef- Since the SHC u s H is proportional to ((3pf) 2 and 
U/ep within our calculation scheme, the vertical axis is 
set to be <j s h ([3pf)~ 2 {£f /U), and is normalized by the 
universal value e/(8n). a s H has a maximum value along 
the line a = 2\ppT = 1 at each A and the maximum 
value depends on \og(A/pp). This conserved spin current 
is nonzero even for the Rashba model in two dimensions. 

The above results give a hint to look for systems show- 
ing the spin Hall current. One important feature is that 
the Hamiltonian should involve s z . Luttinger model^ sat- 
isfies this condition while the cubic Rashba model does 
not. Therefore the complete confinement of the electronic 
motion along one direction is not desirable. We have as- 
sumed that the spin-orbit interaction is unchanged by 
disorder. In reality, the impurity potential v(r) induces 
a spin-orbit coupling as (p x Vv) ■ s. This is also ex- 
pected to contribute to the extrinsic spin Hall effect^ 
This effect is beyond the scope of the present paper. 

In conclusion, we derived a general exact formula of the 
spin Hall conductivity for the two kinds of spin current: 
(a) the product of spin and velocity operators and (b) 
the effective conserved spin current. The conditions for 
the non-zero spin Hall current has been clarified and are 
applied to the Rashba and cubic Rashba models. 
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